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First simultaneous multiwavelength observations of the BL Lac object 1ES1959+650 
in a steady state with MAGIC and Suzaku/Swift 
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Abstract: Simultaneous multiwavelength observations were conducted for the BL Lac object 
1ES 1959+650 in a steady state in May 2006 with the MAGIC telescope and the X-ray satellites Suzaku 
and Swift. Swift can also provide multi-filter photometry in the UV-optical band. The source was clearly 
detected in all observed energy bands, from the optical to TeV. With respect to previous observations the 
source was in a low state in the very high energy (VHE) band (~10% Crab flux above 300 GeV) but in 
a relatively high state in the X-ray band. The light curves showed rather stable activities, with no signif- 
icant variability in the VHE 7-ray emission and small variability (~ 10% amplitude) in the X-ray band. 
The observed spectral energy distribution in the steady state can be described by a one-zone synchrotron 
self-Compton model. 



Introduction 

Multiwavelength observations of the BL Lac ob- 
jects can provide information on the acceleration 
mechanism and the properties of the jet. Since the 
BL Lac objects are highly variable sources, simul- 
taneous multiwavelength observations are essen- 
tial tools to test the various emission models. At 
present, it is unclear whether the observed emis- 
sion is entirely due to accelerated electrons, as in 
synchrotron self-Compton (SSC) models [e.g.] [5], 
or whether the high-energy emission is due to sec- 
ondary electrons produced by accelerated protons 
and ions [e.g.][12]. 

1ES 1959+650 (RA = 19:59:59.95, Dec = +65:08:54 
[J2000], z = 0.047) is categorized as a High- 
frequency BL Lac object (HBL) where the 
synchrotron peak lies in the X-ray regime [14]. 
The first TeV 7-ray signal from this object was 
reported in 1998 [13]. Recent results from the 
observations in 2004 with the MAGIC telescope 
show a signal at the 8 a level and a flux of about 



20% Crab (> 300 GeV) from 6.5 hours of observa- 
tion [2]. This source is particularly interesting for 
multiwavelength observations because a so-called 
"orphan" flare in TeV band without any X-ray 
counterparts was observed in 2002. This activity 
can not be explained by conventional one-zone 
SSC models [10]. 

Observations 
MAGIC Observations 

The MAGIC (Major Atmospheric Gamma Imag- 
ing Cherenkov) telescope is an Imaging Atmo- 
spheric Cherenkov Telescope (IACT) with a 17-m 
diameter dish, located on the Canary Island of La 
Palma (28.2° N, 17.8° W, 2225 ma.s.l.). The tele- 
scope is operating at a 7-ray trigger threshold of 
~ 50 GeV and a spectral threshold of —100 GeV. 
The telescope parameters and performance are de- 
scribed in detail in [3], 



MWL OBSERVATIONS FOR 1ES 1959+650 WITH MAGIC AND SUZAKU/SWIFT 



e 2 plot | 



Suzaku Observations 



« 400 



350 ft 

300 

250 



150 
100 
50 



1ES1 959+650 (z = 0.047) 



e 2 cut 



2006 May 21-27 

Effective observation time: 14.3 h 
zenith angle: 36-43.5 deg 



200 if +!++' 



preliminary 



Background events: 603 +- 14 
Excess events: 324 +- 34 

Excess significance: 10.4 o 



0.05 0.1 0.15 0.2 0.25 0.3 

e 2 [deg 2 ] 



Figure 1: 6 2 distribution of 1ES 1959+650 above 
350 photoelectron (corresponding to the 7-ray 
peak energy of —400 GeV) . A cut of 9 2 < 
0.03 deg 2 (dashed blue line) yields the signal 
of 324 events over 603 normalized background 
events, corresponding to a 10.4 a excess. 



Suzaku [11] is a X-ray satellite with four X-ray 
Imaging Spectrometers (XIS) and a separate Hard 
X-ray Detector (HXD). The XIS are sensitive in 
the 0.2-10 keV band with two types of CCDs com- 
posed of front-illuminated CCDs (for XIS0, 2 and 
3) and a back-illuminated type (XIS1). XIS1 is 
more sensitive below 2 keV. The HXD is a silicon 
PIN diode array and is the most sensitive detector 
in the 10-70 keV band thanks to the good noise 
shielding. 

The Suzaku satellite was pointing 1ES 1959+650 
between 01:13:23 of May 23rd and 04:07:24 of 
May 25th, 2006 (UT). The HXD/PIN light curve 
shows a rapid increase of the noise after about 100 
ksec (possibly due to the in-orbit unexpected radi- 
ation 1 ) and the data after this event cannot be used 
for the analysis. After screening the total net expo- 
sure time for XIS and HXD are 99.3 ksec and 40.2 
ksec, respectively. According to the RXTE/ASM 
light curve of May 2006 2 , the source was in a rel- 
atively high state during the Suzaku pointing. 



1ES 1959+650 was observed with the MAGIC tele- 
scope from 21st to 27th of May, 2006. The zenith 
angle during the observations was in the range 
from 36° to 43.5°. Observations were performed 
in wobble mode [6] with the object offset by 0.4° 
from the camera center. After the quality selection 
of the data the total effective observation time was 
14.3 hours. The analysis was performed using the 
standard MAGIC analysis software. Details of the 
analysis chain can be seen in [3]. In order to de- 
termine the significance of the 7-ray signal and the 
number of excess events, a final cut with 8 2 pa- 
rameter was used. The is the angular distance 
between the source position in the sky and the re- 
constructed arrival position of the air shower which 
is estimated using "DISP" method [7]. A cut of 
2 < 0.03 deg 2 yields an excess of 324 events over 
603 normalized background events, corresponding 
to a 10.4 a excess as shown in figure 1 . Only events 
> 350 photoelectrons were used, corresponding to 
the peak of the energy distribution at about 400 
GeV for 7-ray. 



Swift Observations 

The Swift satellite [8], with its easy and flexible 
scheduling, can be optimally employed for the ob- 
servation of bright blazars. It carries three instru- 
ments. The Burst Alert Telescope (BAT) is opti- 
mized for 15-150 keV with high sensitivity while 
the X-Ray Telescope (XRT) is sensitive in the 0.3 - 
10 keV band. In addition a UV/optical telescope 
can provide data in the 170-600 nm band. The 
source was also observed with Swift around the 
Suzaku pointing 19th, 21st and from 23rd to 29th 
of May. 

Results and Discussions 
Light curves 

The light curves in VHE 7-ray and X-ray are 
shown in figure 2. The average integrated flux 
above 300 GeV is (1.27 ±0.16) x 10 _n cm- 2 s _1 
(x 2 /d.o.f = 8.5/6), which corresponds to about 10% 



1. see http://www.astro.isas.ac.jp/suzaku/log/hxd/ 

2. see http://heasarc.gsfc.nasa.gov/xte_weather/ 
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of the Crab Nebula flux [3]. No significant strong 
variability can be seen in VHE 7-ray emission. 
In the X-ray band a flare of small amplitude (~ 
10%) with a rising time of T t ~ 10 5 s could be 
observed. The hardness ratio shows rather small 
spectral variation, except for a the sudden drop at 
t ~ 1.5 x 10 5 s. 

Spectra 

Since the source was stable both in the VHE 7- 
ray and the X-ray bands we extracted the spectra 
using the whole dataset during the campaign for 
each band. 

Figure 3 shows the measured spectrum by 
the MAGIC telescope. It is well described 
by a simple power law from 150 GeV to 
3 TeV with dN/dE = (2.7 ± 0.3) x 
10- 12 {E/l TcV )- 2 - 58±0 - 18 [TeV- 1 s- 1 cm- 2 ]. 
Compared to the previous MAGIC measurement 
of 1ES 1959+650 in a steady state in 2004 [2], the 
observed flux in 2006 was 60% of the 2004 result 
while the photon indices are comparable. 

In X-ray, XIS spectra were extracted for t < 10 5 
s in order to perform a joint XIS and HXD/PIN 
fit. The model with three power-laws is re- 
quired to obtain reasonable fitting results. The 
best fit parameters with associated errors are: 

r x = 1.94 ± 0.01, r 2 = 2.195 ± 0.003, r 3 = 

2.75^5, £br,l = 1.83±0.03,£ b r.2 = 16^ with 

X 2 /d.o.f = 4065.0/4139. 

Spectral Energy Distributions 

In figure 4 we report the spectral energy distribu- 
tion (SED) of 1ES 1959+650 as measured at the 
end of May 2006, together with other historical 
data [see the figure caption for the detailed infor- 
mation]. The absorption of VHE 7-ray emission 
by the extra-galactic background light (EBL) are 
corrected using "Low model" of [9]. 
Assuming the uniform injection of the electrons 
throughout a homogeneous emission region we 
model the observed SED by using a one-zone SSC 
model developed by [16, 17]. Briefly, a spheri- 
cal shape (blob) is adopted for the emission region 
with a radius R, filled with a tangled magnetic field 
with intensity B . The emission region is moving 
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Figure 2: Light curves of 1ES 1959+650 during 
this campaign. [Upper]: X-ray light curves ob- 
tained by the Suzaku XIS 1 detector. Top to bottom: 
Count rates in the 0.2-2 keV band, in the 2-10 keV 
band and the hardness ratio (2-10 keV/0.2-1 keV). 
[Lower]: Diurnal light curve of VHE 7-ray (> 300 
GeV) measured by the MAGIC telescope. A dotted 
horizontal line represents the average flux. Dotted 
vertical lines show the observation window of the 
Suzaku pointing. 
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1ES1 959+650 measured spectrum [21-27 May 2006] 
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Figure 3: Differential energy distribution of 
1ES 1959+650 averaged over the whole dataset of 
the campaign in 2006. The blue solid line repre- 
sents a power-law fit to the measured spectrum. 
The fit parameters are listed in the figure. For com- 
parison, the measured MAGIC Crab spectrum [3] 
is shown by a red dashed line. 
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Figure 4: Overall SED of 1ES1959+650. Optical- 
UV data are from REM (cyan triangle) and 
UVOT/Swift (blue triangles). We report the av- 
erage Suzaku spectrum (red) and the swift spec- 
tra taken on May 24 [higher] and May 29 [lower] . 
Green points (filled circles) report the observed 
MAGIC spectrum, while the red points (empty 
triangles) have been corrected for the absorption. 
Historical data are taken from [14] (radio-optical- 
X-rays), [4, 10] (X-rays), and [1] (TeV). The line 
refers to the best fit SSC model for the campaign 
data in 2006 using the code developed by [16, 17]. 
See text for the fit parameters. 

at a relativistic speed (3c with a bulk Lorentz fac- 
tor r and the viewing angle of 6. The beaming 
(Doppler) factor is given by 6 = 1/[T(1— /3cos#)], 
with 5 ~ T for 9 ~ 1 /Y. The electron distribu- 
tion is described by a smoothed broken power-law 
JV( 7 ) = if 7 - Sl (l + 7/7b) Sl " S2 exp(- 7 / 7max ) 
where if is a normalization factor, 7 b is the break 
Lorentz factor of electrons, and s% and s 2 are the 
spectral indices below and above the break, respec- 
tively. The best fit SSC model was achieved using 
the following parameters: S = 35, R — 5 x 10 15 
cm, B = 0.12 G, K = 2 x 10 3 cm" 3 , 7min = 1, 
7b = 8 x 10 4 , 7max = 7 x 10 5 , si = 2 and 
s 2 = 3.5. 

Summary 

For the first time, in May 2006, we successfully 
performed simultaneous multiwavelength observa- 
tions in the VHE 7 -ray, X-ray and optical-UV 



bands for 1ES 1959+650 in a quiescent state. The 
source was clearly detected in all observed en- 
ergy bands. The results indicate that the activity 
of the source was rather stable. The VHE 7 -ray 
flux showed a low state of activity with about 10% 
Crab flux unit above 300 GeV and a photon index 
of -2.58 ± 0.18. In the X-ray band the flux was 
relatively high and only small variability (~ 10%) 
can be seen. The SED with the simultaneous spec- 
tral data in the steady state can be well described 
by the one-zone SSC model. More details will be 
available in [15]. 
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